LATERAL TUNNELING THROUGH THE CONTROLLED BARRIER BETWEEN 
EDGE CHANNELS IN A TWO-DIMENSIONAL ELECTRON SYSTEM 



A.A. Shashkin, V.T. Dolgopolov, E.V. Deviatov, B. Inner f, A.G.C. Haubrich f, J.P. Kotthaus f, M. Bichler t, 

W. Wegscheider £ 

Institute of Solid State Physics, Chernogolovka, Moscow District 142432, Russia 
\Ludwig-Maximilians-Universitat, Geschwister-Scholl-Platz 1, D-80539 Miinchen, Germany 
% Walter Schottky Institut, Technische Universitat Miinchen, D-85748 Garching, Germany 



We study the lateral tunneling through the gate- voltage-controlled barrier, which arises as a result 
of partial elimination of the donor layer of a heterostructure along a fine strip using an atomic 
force microscope, between edge channels at the depletion-induced edges of a gated two-dimensional 
electron system. For a sufficiently high barrier a typical current-voltage characteristic is found to be 
strongly asymmetric and include, apart from a positive tunneling branch, the negative branch that 
corresponds to the current overflowing the barrier. We establish the barrier height depends linearly 
on both gate voltage and magnetic field and we describe the data in terms of electron tunneling 
between the outermost edge channels. 



Recently there has arisen much interest in the lateral 
tunneling into the edge of a two-dimensional electron sys- 
tem (2DES), which is related not only to the problem of 
integer and fractional edge states in the 2DES but also to 
that of resonant tunneling and Coulomb-blockade [Q-Q. 
The tunneling regime was identified by exponential de- 
pendences of the measured current on either source-drain 
voltage Jl|-§1 or magnetic field || . For producing a tunnel 
barrier a number of methods were used: (i) gate voltage 
depletion of a narrow region inside the 2DES 0^J^1 ; (ii) 
focused- ion-beam insulation writing ||] ; (iii) cleaved-edge 
overgrowth technique As long as the tunnel barrier 
parameters arc not well controllable values, it is impor- 
tant that using the first method one can tune the barrier 
on the same sample. In contrast to vertical tunneling 
into the bulk of the 2DES at a quantizing magnetic field, 
when the 2DES spectrum shows up |^|,^), in the lateral 
tunneling electrons can always tunnel into the Landau 
levels that bend up at the edge to form edge channels 
where they intersect the Fermi level, i.e., the spectrum 
gaps are not seen directly in lateral tunneling. Instead, it 
reflects the edge channel structure and density of states. 
For both the integer and fractional quantum Hall effect, 
a power-law behaviour of the density of states at the 
2DES edge is expected. This can interfere with the bar- 
rier distortion at an electric field in the nonlinear regime 
of response and, therefore, results of lateral tunneling 
experiments obtained from the measurements of current- 
voltage curves || should be treated with care. 

Here, we investigate the lateral tunneling in the nar- 
row constrictions in which, along a thin strip across, 
the donor layer of a GaAs/AlGaAs heterostructure is 
partly removed using an atomic force microscope (AFM) . 
A controlled tunnel barrier is created by gate depletion 
of the whole of the sample. The well-developed tunnel- 
ing regime is indicated by strongly asymmetric diode-like 
current-voltage characteristics of the constriction which 



are sensitive to both gate voltage V g and normal magnetic 
field B. The behaviour of the tunneling part of current- 
voltage curves points to electron tunneling between the 
outermost edge channels. 

The samples are triangular constrictions of a 2D elec- 
tron layer with different widths W — 0.7, 0.4, 0.3, and 
0.2 /im of the thinnest part, see Fig. |l|(a). These are made 
using a standard optical and electron beam lithography 
from a wafer of GaAs/AlGaAs heterostructure with low- 
temperature mobility \i = 1.6 x 10 6 cm 2 /Vs and carrier 
density n s = 4 x 10 11 cm~ 2 . Within each constriction 
the donor layer is removed along a fine line by locally 
oxidizing the heterostructure using AFM induced oxida- 
tion jl(| . This technique allows one to define 140 A wide 
oxide lines of sufficient depth and oxide quality so as to 
partly remove the donor layer and, therefore, locally de- 
crease the original electron density. The whole structure 
is covered with a metallic gate, which enables us to tune 
the carrier density everywhere in the sample. When de- 
pleting the 2D layer the oxidized regions get depopulated 
first, resulting in the creation of tunnel barriers. Poten- 
tial probes are made to the sample to allow transport 
measurements. 

For the measurements we apply a dc voltage, V s d, be- 
tween source (grounded) and drain contacts of one of 
the constrictions modulated with small ac voltage with 
amplitude V ac = 40 fiV and frequency / = 20 Hz. A 
gate voltage is applied between the source and the gate. 
We measure the real part of the ac current, which is 
proportional to the differential conductance dl/dV, as a 
function of bias voltage V s d (I — V characteristics) us- 
ing a home-made I — V converter and a standard lock-in 
technique. The behaviour of I ~ V characteristics is in- 
vestigated with changing both gate voltage and magnetic 
field. The measurements are performed at a temperature 
of about 30 mK at magnetic fields of up to 14 T. The re- 
sults obtained on different constrictions are qualitatively 
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FIG. 1. (a) Top view on the sample (top); and a blow-uf 
of one of the constrictions upon etching the oxidized part o: 
the mesa as performed solely for visualization purpose (bot- 
tom), (b) Gate voltage dependences of the electron density ir 
the oxidized (squares) and unoxidized (circles) regions of the 
2DES. An example of the magneto-conductance in the barriei 
region is shown in the inset. The value of n a is extracted fron 
the slope of dashed lines with 10% uncertainty. 

similar. 

To characterize the sample we extract the gate voltage 
dependence of the electron density from the behavioui 
of magneto-conductance plateaux in the barrier regior 
and in the rest of the 2DES (Fig. |(b)). The analysis 
is made at high fields where the size-quantization-causec 
effect of conductance plateaux in narrow constrictions i; 
dominated by magnetic field quantization effects jllj . As 
seen from Fig. ^(b) , if the barrier region is depopulatcc 
(V g < Vth), the electron density in surrounding areas is 
still high to provide good conduction. The slopes of the 
dependences n s (V g ) in the oxidized region and in the res1 
of the 2DES turn out to be equal within our accuracy 
The distance between the gate and the 2DES is deter- 
mined to be d ~ 570 A; as the corresponding growth 
parameter is about 400 A, the 2D layer thickness con- 
tributes appreciably to the distance d. We have found 
that even in the unoxidized region the electron density 
at V g — can be different after different coolings of the 
sample because of insignificant threshold shifts: it falls 
within the range 2.5 x 10 11 to 4 x 10 11 cm -2 and is always 
higher compared to the barrier region. 

A typical I — V characteristic of the constriction in the 
well-developed tunneling regime is strongly asymmetric 
and includes an overflowing branch at V s d < and the 
tunneling branch at V s d > 0, see Fig. ||(a). The tunnel- 
ing branch is much smaller and saturates rapidly in zero 
B with increasing bias voltage. The onset voltages Vq 
and Vt for these branches are defined in a standard way 
as shown in the figure. The tunneling regime can be at- 
tained by both decreasing the gate voltage and increasing 
the magnetic field, as is evident from Fig.g(a). We check 



that the shape of / — V characteristics is not influenced 

by interchanging source and drain contacts. Hence, the 
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FIG. 2. (a) I— V curves at different gate voltages and mag- 
netic fields. The cases of B = and B ^ correspond to two 
coolings of the sample as compared in Fig. ^|(a). W = 0.4 /im. 
(b) A sketch of the 2D band bottom in the barrier region for 
different source-drain biases V s d- 

To understand the asymmetry origin let us consider 
a gated 2DES containing a potential barrier of approxi- 
mately rectangular shape with width L 3> d in zero mag- 
netic field. The 2D band bottom in the barrier region co- 
incides with the Fermi level Ep of the 2DES at V g equal 
to the threshold voltage Vth* Since in the barrier region 
for V g < Vth an incremental electric field is not screened, 
the 2D band bottom follows the gate potential so that the 
barrier height is equal to —eAV = e(Vth — V g ), where 
— e is an electron charge (Fig. 0(b)). Applying a bias 
voltage V s d leads to shifting the Fermi level in the drain 
contact by —eV s d- Because of gate screening the voltage 
V s d drops on the scale of the order of d near the boundary 
between barrier and drain and so the barrier height on 
the source side does not practically change, see Fig. ||(b). 
If V s d reaches the onset voltage Vo — AV g , the barrier 
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on the drain side vanishes and electrons start to overflow 
from drain into source. In contrast, for V s d > the elec- 
tron tunneling through the barrier from source into drain 
is only possible. With increasing V s d above —AV g , the 
tunneling distance diminishes and the barrier shape be- 
comes close to triangular. Within the triangular barrier 
approximation, in the quasiclassical limit of small tunnel- 
ing probabilities, it is easy to deduce that the derivative 
of the tunneling current over bias voltage is expressed by 
the relation 
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where <tq k —(e 2 /h)AV g W/V s d^F, m = 0.067mo (mo is 
the free electron mass), and Xp is the Fermi wave-length 
in the source. Obviously, the tunneling current is domi- 
nated by electrons in the vicinity of the Fermi level, and 
the tunneling distance L T = —AV g L/V s d should satisfy 
the inequality d <C Lt < L. In accordance with Eq. ([!]), 
the expected dependence of the tunneling onset voltage 
Vt on gate voltage is given by Vt oc (— AV g ) 3 ^ 2 . 

As seen from Fig. [|(a) , the expected behaviour of both 
Vo and Vt with changing V g is indeed the case. The 

dependences Vo(V g ) and V^ 3 (V g ) are both linear; the 
slope of the former is very close to one. Extensions of 
these straight lines intercept the V^-axis at slightly dif- 
ferent voltages, which points out that the triangular bar- 
rier approximation is good. The threshold voltage Vth for 
the 2DES's generation in the barrier region, which is de- 
fined as a point of vanishing Vo (Fig. ||(a)), is coincident, 
within experimental uncertainty, with the value of Vth 
determined from the analysis of magneto-conductance 
plateaux (Fig. |l| 




FIG. 3. (a) Change of the onset voltages Vo and Vt as de- 
fined in Fig. |^(a) with V g at B — 0; and (b) the fit (dashed 
lines) of I — V curves (solid lines) by Eq. (jl|) with the pa- 
rameters L = 0.6 /im, (j = 38 Mfi -1 , Vth = -0.4 mV; 
W — 0.4 /jm. In case (a), the data marked by filled triangles 
are obtained for the same cooling of the sample as the ones at 
B — in Fig. ^(a) and in case (b) , whereas the open triangles 
correspond to the B / data of Figs. §(a) and ^[as measured 
for the other cooling. 



Fitting the set of I - V characteristics at different V„ v a > v th, where there is no tunnel barrier in zero B, 



by Eq. (0) with parameters L, Vth, and Co is depicted in 
Fig. ||(b). The dependence of <7o on AV g and V s d is ig- 
nored on the background of the strong exponential depen- 
dence of dl/dV. Although three parameters are varied, 
the fit is very sensitive, except for ctq, to their variation 
because of the exponential behaviour of I — V character- 
istics. One can see from Fig. |](b) that the above model 
describes well the experiment at zero magnetic field. As 
expected, the determined parameter L = 0.6 fim is much 
larger than d, i.e. the barrier shape at V s d = is ap- 
proximately rectangular, and the value of Vth is close to 
the point where Vo (and Vt) tends to zero (Fig. 0(a)). 
The similar results are obtained on the other two con- 
strictions. Besides, we find that the coefficient (To for 
different constrictions does not scale with the constric- 
tion width W. This probably implies that the tunnel 
barriers even with submicron lengths are still inhomoge- 
neous, which, however, does not seem crucial for the case 
of exponential I — V dependences. 

Having tested that we deal with the controlled tunnel 
barrier we investigate the tunneling at a normal magnetic 
field that gives rise to emerging the tunnel barrier in a 
similar way to gate depletion (Fig. ||(a)). At constant 



the magneto-conductance a obeys the 1/B law at weak 
fields and drops exponentially with B in the high-field 
limit, signaling the tunneling regime. Fig. |J(a) presents 
the magnetic field dependence of the onset voltage Vo 
that determines the barrier height. As seen from the 
figure, the change of the barrier height —eVo with B is 
very close to hui c /2, which points to a shift of the 2D 
band bottom by half of the cyclotron energy. 

For describing the tunneling branch of I — V character- 
istics we calculate the tunneling probability in the pres- 
ence of a magnetic field. This is not so trivial as at B = 
because electrons tunnel through the magnetic parabola 
between edge channels at the induced edges of the 2DES. 
In the triangular barrier approximation one has to solve 
the Schrodinger equation with the barrier potential 



U{x) = -^{x 



eV sd -~eAV g , < x < L, (2) 



where lo c is the cyclotron frequency, I is the magnetic 
length, and eV s d is larger than the barrier height in the 
magnetic field. An electron at the Fermi level in the 
source tunnels through U(x) from the origin to a state 
with orbit centre Xq such that < xq < L. If the 
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barrier potential is dominated by the magnetic parabola 
(i.e., the magnetic length is the shortest), the problem 
reduces to the known one to find energy levels in the 
shifted parabolic potential as caused by the linear term 
in Eq. (0) . The value of xq is determined from the con- 
dition of coincidence of a Landau level in the potential 
U(x) with the Fermi level in the source. If the lowest 
Landau level is regarded only and the spin splitting is 
ignored, we get the minimum tunneling distance to the 
outermost edge channel in the drain 
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The first term in brackets in Eq. (|j|), which is dominant, 
is large compared to one. Knowing the wave function 
of the lowest Landau level in the potential U (x) and ne- 
glecting the last term in Eq. (||) we obtain for the shape of 
I — V characteristics near the onset, where the tunneling 
probability is small, 
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Here <7b is the pre-factor which can be tentatively ex- 
pected to be of the same order of magnitude as oq. From 
Eq. (^) it follows that at sufficiently strong magnetic 
fields the tunneling onset voltage Vr is related to the 
barrier height as Vrl oc Huj c /2 — eAV g , which is consis- 
tent with the experiment (Fig. ^(a)). The solution (^) 
includes the case eAV g > when a tunnel barrier is ab- 
sent at zero magnetic field but arises with increasing B. 
This occurs apparently because of depopulation of the 
barrier region in the extreme quantum limit of magnetic 
field. 

Fig. ||(b) displays the fit of / — V characteristics at 
different magnetic fields by Eq. ([|) with parameters L, 
Vth, and ub- The optimum values of L = 0.6 /Ltm and 
Vth — —1.5 mV are found to be very close to the ones 
for the B = case as determined for the same range of 
barrier heights, see Fig. ||(b). Although this fact sup- 
ports our considerations, these are not quite rigorous to 
discuss the considerable discrepancy between the fore- 
exponential factors with and without magnetic field. 

The observed behaviour of / — V characteristics with 
magnetic field in the transient region where their asym- 
metry is not yet strong (Fig. |^(a)) is similar to that of 
Refs. |,|. Over this region, which is next to the scope of 
the exponential I — V dependences at higher magnetic- 
field-induced tunnel barriers, our I — V curves are close 
to power-law dependences, as was discussed in Ref. [||. 
There is little doubt that it is very difficult to analyze and 
interpret such I — V curves without solving the tunneling 
problem strictly. We note that the peak structures on the 
tunneling branch oi I — V characteristics (see Figs. ||(a) 
and ||(b)) persist at relatively low magnetic fields and are 
very similar to those studied in Ref. M . These may be a 
hint to resonant tunneling through impurity states below 
the 2D band bottom. 
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FIG. 4. (a) Behaviour of the onset voltages Vo and Vr with 
magnetic field; and (b) the fit (dashed lines) of I — V curves 
(solid lines) by Eq. (0) with the parameters L = 0.6 pm, 
cr = 1.3 MSTT 1 , Vth = -1.4 mV; W = 0.4 /im, V a = 0. 
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